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ABSTRACT 

System theorists understand that the same mathematical objects 
which determine controllability for nonlinear control systems of 
ordinary differential equations also determine hypoel 1 iptici ty for 
linear partial differential equations. Moreover, almost any study of 
o.d.e. systems begins with linear systems. It is remarkable that 
Hormander’s paper on hypoellipticity of second order linear p.d.e.’s 
starts with equations due to Kolmogorov, which we show are analogous to 
the linear o.d.e. * s. Eigenvalue placement by state feedback for a 
controllable linear system can be paralleled for a Kolmogorov equation 
if an appropriate type of feedback is introduced. Results concerning 
transformations of nonlinear systems to linear systems are similar to 
results for transforming a linear p.d.e. to a Kolmogorov equation. 
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I INTRODUCTION 


We consider second order linear partial differential equations of 
the forms. 
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where the coefficients are <€ in some open set in IR containing the 
origin. Here we assume the matrix (A^(x)) is symmetric, is positive 
semidef ini te , and has constant rank m . If we consider f as the 
input to the p.d.e and u as its output, then we can replace f by 
f - C(x)u, so we consider C(x) = 0. 

Of particular interest to us is the spatial partial differential 
equation 


(3) 


2 A (x) 
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du 
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f(x). 


with the same assumptions on 
such an equation as 


(Aj k ( x ) ) 


as above. 


Hormander [1] writes 


m 2 2 

(4) 2 Xu + Xou = f(x), 

j=l J 


where Xo. Xi, . . . ,X are 


m 


00 

*€ vector fields on IR f and Xi,X 2 ..-*.X 
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are assumed to be linearly independent. Equations (3) and (4) are said 
to be hypoelliptic if f being < €°° implies that u is <6 . 

Results concerning hypoellipticity of (4) can be compared with 
corresponding results concerning controllability of the nonlinear system 

m 

(5) x = Xo(x) + 2 X.(x)v., 

j=l J J 

where vi,V 2 v are the controls. The linear version of (5) (in 

m 

which the usual roles of A and B are reversed to fit the standard 
upcoming p.d.e. notation) is 

(6) x = Bx + Av, 

where A and B are appropriate matrices. 

Our purpose in this paper is to introduce a particular linear 
partial differential equation (called a Kolmogorov equation) that 
relates to the p.d.e. (3) as the linear system (6) relates to the o.cl e 
(5). For this particular p.d.e. we shall introduce an appropriate 
feedback that allows "eigenvalue placement" if the equation is 
hypoelliptic. We study the effect of this feedback on the spatial 
Fourier transform of the solution u. We also mention the problem of 
transforming (by state coordinate changes and feedback) the linear 
p.d.e. (3) to a Kolmogorow equation as one would transform the nonlinear 
system (5) to a controllable linear system (6). 

We remark that the results of Hormander on hypoellipticity have 
been applied in systems theory by Elliott [2], [3]. 



II. P.D.E. *s AND FEEDBACK 


An equation of special interest is 


( 7 ) 


a 2 u 
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;k=f j7k=i J k j ax k 


= f. 


where A = (a., ) and B = (b ) are constant matrices, A is 

Jk 

symmetric, positive semidef ini te of rank m, and the matrix 
(A,BA.B 2 A, . . . ,B n *A) has rank n (this last condition implies 
hypoellipticity) . Such a p.d.e. is called a Kolmogorov equation . In 
fact if - is added as in (1), Kolmogorov [4] indicates that certain 
probability density functions satisfy such equations. 

By using the input f we can introduce an appropriate type of 
feedback for the equation (7), This feedback should leave invariant the 
principal symbol of (7) and the hypoel 1 iptici ty condition. Writing (7) 
in Hormander’s vector field notation yields 


m 2 

(8) 2 X. u + Xo u = f, 

j =1 3 


where Xi,X 2 X^ are linearly independent and constant coefficient 

and Xo is a linear vector field. Our 1 i near f eedback takes the form 


m _ 

(9) f = 2 k .x X .u, 

j=l 3 3 


where each k^ is a lxn matrix of constants. We can feed back a sum 
of terms involving a linear combination of x variables times a vector 
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field Xj applied to the solution u. Each Xj must be a vector field 
whose square appears in (8). 

2 

Example : Consider the p.d.e. on 1R 
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This is of the form (7) with 
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and of the form (8) with m=l and 
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where ao and at are constants, then (10) becomes 


(13) 


5 u 5u , , 5 ■=■ 

+ X 2 - (a 0 xi + a 1 X 2 ) = f. 

5xf 5xi 5x2 


Moreover , 


A remains unchanged and B moves to 
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Neither the principal symbol or hypoelliptici ty of (10) is changed by 
the linear feedback. We can choose ao and ai to yield any desired 
characteristic polynomial for B. 

We want to consider the effect of feedback of the form (9) on the 
p.d.e. (7). Following the argument in [1] we take Fourier transforms 

A 

with f = (fi,f 2 f n ) being the transform variables. Then (with u 

being the Fourier transform of u) we find 


n 


(15) 


- A( f ,nu-^ iVk |- 


= f. 


This is a first order linear p.d.e. which can be solved by the method of 
characteristics for a given noncharacteristic set of initial conditions. 
The characteristic curves are determined by (B ' denotes B transpose) 


(16) £ = B 'f 


A 

and u must satisfy 


(17) 57 = A(f.f) - f. 


Here r denotes the parameter along the characteristic curves. 

Hence the linear feedback in which the B matrix is altered simply 

/v 

changes the characteristic curves used in solving for u. Equation (17) 
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is actually unaltered. 


Given the Kolmogorov equation (7) with matrices A and B, we can 
form the controllable linear control system (6) 

x = Bx + Av. 


We define the Kolmogorov indices £1 > £2 > . . . 2 of (7) to be the 

Kronecker indices ki > *2 > . . . > k of this linear system. Canonical 

“ m 

forms for (7) which parallel canonical forms for linear systems can also 
be derived. 

We state our main result without proof. 


Theorem Assume that the matrices 
p.d.e. (7) 
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satisfy rank [A, BA B A] = n. Then the eigenvalues of the B 

matrix can be arbitrarily placed (with complex eigenvalues occuring in 
conjugate pairs) by linear feedback. 

If rank A = m is n, then equation (7) is elliptic. Then linear 
feedback can be used to eliminate all first order terms. This compares 
to x = Bx + Av with a control for each state. 

Our Kolmogorov equation (7) has constant coefficient second order 
part and linear varying first order part. However, our original partial 
differential equation (3) has variable coefficients. An interesting 
problem is to derive necessary and sufficient conditions to transform 
(in a designated sense) the general hypoelliptic equation (3) to a 
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Kolmogorov equation (7). This should compare to results for moving a 
nonlinear control system of o.d.e’s to a controllable linear system [5], 

[6]. [7], [8]. [9]. [10]. 

Of course, partial differential equations of order higher than two 
can be considered. Also the implementations of the results of this 
paper by finite difference and finite element schemes when only point 
sensors and actuators are involved is an open problem. 


Ill CONCLUSION 

We have drawn parallels between control theory for linear o.d.e’s 
and the Kolmogorov p.d.e.’s. An appropriate type of linear feedback for 
p.d.e.’s was introduced. A transformation theory for more general 
partial differential equations is presently being researched. 
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